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^ ' A non-linear Boltzmann equation|3 describing the time evolution of a partonic sypi- 

O ' tern in the central rapidity region after a heavy ion collision is solved numericallycl. 

' A particn^ model of the coUinear logarithmic divergences due to small angle 
scatteringau is employed in the numerical solution. The system is followed until it 

y—^ , reaches kinetic equilibrium where the equilibration time, temperature and chemical 

1 potential are determined for both RHIC and LHC. 

> : 

^-H , 1 Introduction 

O ; 

\ Heavy ion collisions are already taking place at the Relativistic Heavy Ion 

■ Collider (RHIC) and are planned to take place in a few years at the Large 
(«~^ \ Hadron Collider (LHC). The question remains whether the hot and dense 

' matter formed after the collision equilibrates to form the quark gluon plasma. 

\ After a heavy ion collision, an initial distribution of partons, way out of 

^ • equilibriumEl, is formed and strongly influences the subsequent interactions 

I \ of the partons with each other which drive the partonic system towards equi- 

librium. Whether or not the system equilibrates depends on if the time it 
(— I I takes to reach equilibrium is shorter than the hydrodynamic expansion time. 

' The small x Fock states in nuclei responsible for multi-particle produc- 

. \ tion at central rapidities have distributions that are described by a classical 

■ effective field theory (EFTlEl. The classical distributions in a single nucleus 
H \ can be solved analyticallyO'S. The classical gluon distribution falls off as 1/ /c^ 

at large transverse momentum kt but saturates at smaller kt] growing loga- 
rithmically with kt. Non-perturbative, numerical solutions of the Yang-Mills 
equations have determined exactly initial number and energy distribution of 
gluons after a coUisioncI and provide an initial condition for the single parti- 
cle gluons distributions in a transport ecmation. For simplicity, an idealized 
approximation for the gluon multiplicity^ is used. Long ago the parton evo- 
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lution was described by a Boltzmann equation under the assumption of boost 
invariance and solved in the relaxation time approximationEI. RecenJjy, a 
non-linear Boltzmaiin description of a partonic system, was formulatedcl and 
solved numericallya, with no approximation, for a particular pmodel of the 
collinear logarithmic divergences due to small angle scatteringErQ. The time 
and temperature for the system to reach kinetic equilibrium for both RHIC 
and LHC were calculated. This Boltzmann description takes into account 
number conserving 2 — > 2 processes. 2 — > 3 processes, which may have a 
large effectlliJ, were not considered. 

The distribution of gluons at the very early times after a high energy 
heavy ion collision is described by the bulk scale Qs of gluon saturation in the 
nuclear wavefunction. The scale Qs is the only scale in the problem, and all 
observables calculated depend sensitively on it. The gluon system is followed 
until it reaches kinetic equilibrium where the equilibration time, temperature 
and chemical potential are determined In this paper we take as = 0.3 and 
Nc = 3. 



2 The Boltzmann equation for gluons 

The Boltzmann equation is: 



^+«-V/ = C(/), (1) 

where C(/) is the collision integr|al and / = f{t, x,p) is the single particle 
gluon distribution. Following Baym □ we assume that the transverse dimen- 
sion of the collision volume for central collision is large enough for / to depend 
spatially only on the coordinate along the collision axis, z. Next we assume 
that the central rapidity region is invariant under Lorenz± boosts Ii3. Under 
these assumptions we find that v ■ V/ = I'zfl = ~^^r " where / = f{t,p) 
since we can always choose z = due to the assumption of boost invariance. 
The l.h.s of eq. (0) changes accordingly. The collision integral C(/), was 
determined by MuellcrEI. The Boltzmann equation becomes Ela: 



— dt 1 dp — " Mt)L\7p if) + 2\n-i{t)LVp ■ (vf) , (2) 

where, 

f d^P f d?P f 
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where A = 27rag JILi i ^ — p/IpIi 9g — 2(A''^ — 1), n is the number density 
of gluons. Vp' refers to differentiation with respect to p. Integration of (||) 
with respect to (Pp, leads us to the condition tn ~ constant. Furthermore, 
L represents_the logarithmic collinear divergence arising from small angle 
scattering 



where ^) is the average energy and mo is the Debye screening mass 
defined as 



asNc f (Pp 



(5) 



3 Initial conditions 



We solve eq. (||), using the numerical method discussed in a previous paperB. 
The initial gluon multiplicity, dN / (Pp± , for high energy nuclear collisions can 
be calculated along with the formation time of the gluons that are produced 
in the collision. Therefore the initial gluon single particle distributions f{t,p) 
required to solve eq. (|^) are determined: 

(27r)3 dN S{p,) 

f[t,P) = -72 7 ' (6) 

go d^pi_ to 

where is 1.35 GeV"^ for RHIC and 0.65 GeV"^ for LHC. At RHIC 
Qs ~ 1 GeV and at LHC Qs ^ 2 — 3 GeV. Although the initial gluon single 
particle distributions, Eq. (^, were recently determined by Krasnitz and 
VenugonalanEl for simplicity we approximate the p± distribution with a step 
functiontl: 



as-l^cto 

4 Bulk quantities at equilibrium 

Having numerically solved eq. for / as a function of t and p it is useful to 
define bulk quantities of the gluon system, such as the local energy density e 
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and the longitudinal pressure Pl, 

<t)^gol^pfit,p) , P,it)^9al^ff{t,P). (8) 

We note also that the transverse pressure Pt = ^ (e — Pf, ) . Furthermore, 
for Boltzmann statistics, the definition of the entropy density per particle is: 

s^-^^l f{t,p)logf{t,p). (9) 

At kinetic equilibrium 3Pl — SPt — £• Therefore from Bjorken hydro- 
dynamics 113 the energy density e ^ t^^^^ . At kinetic equilibrium, the single 
particle distribution is f{t,p) — exp | p |)) where f3 ~ 1/T. Some useful 

quantities at kinetic equilibrium are: 

n{t) = — 5- exp-, n_i(t) = — -^exp- , e = 3 — 5- exp-. (10) 

From eq.'s (pO|), we note that at kinetic equilibrium (since the number 
density n ^ l/t) T ~ i^^/^, and that /i/T ^ constant with respect to time. 
Furthermore, the temperature at kinetic equilibrium is proportional to the 
energy density per particle, T e/n. 

5 Numerical results 

Fig. 1^ shows, for a, and Qs at RHIC energies, the energy density per particle, 
e/n, times t^/'^ ^ the entropy density per particle and logarithm L versus time. 
As shown, tt^l"^ jn increases monotonically and asymptotically approaches a 
constant, as expected at kinetic equilibrium. The entropy density per par- 
ticle s also monotonically increases and approaches a constant. This is also 
expected since at equilibrium, s = 3 — /^/T. Lastly, fig. ^ shows L (defined 
in eq. ^) as a function of time. Since (p) = e oc T at kinetic equilibrium and 
m|j oc T^, the argument of the logarithm in cq. (^) approaches a constant at 
equilibrium. Fig. ^ shows that i is a constant at equilibrium. We define the 
time where kinetic equilibrium begins to set in as the time igq, where s and 
t^^^t/n reach 10% of their asymptotic value. Having determined the time, iegj 
at equilibrium we determine the the temperature Teg, from ( pO| ) and finally 
the chemical potential ^. 

In this way, for as = 0.3 and Nc — 3, we find the equilibration time at 
RHIC {Qs = 1 GeV), t^q = 3.24 fm, the temperature T^g = 174.27 MeV, 
and the chemical potential ^ = 157.86 MeV. At LHC (Q^ = 2 - 3 GeV), 
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« = 0.3, c = 1 .3, Q. = 1 .0, Nc=3 
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Figure 1. t^/^ times the energy density per particle, the energy density e, and the lon- 
gitudinal and transverse pressures, and Pt, the entropy density per particle s and L 
plotted as functions of time. 



teq = 1.42 ~ 2.36 fm, the temperature Te, = 320.72 ~ 471.69 MeV, and 
H = 249.61 ~ 457.72 MeV. 



6 Conclusion 

In this paper, the Boltzmann equation for gluons was solved numerically for 
the single particle distributions f{t,p). The time evolution of f{t,p) was 
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followed to kinetic equilibrium. The energy density e and entropy density 
per particle s were determined. The results are expected, to change if Bose 
enhancements are included in the Boltzmann equationOiUand a realistic ini- 
tial single paiticle distribution, supplied by the calculation by Krasnitz and 
VenugopalanEl is used as the initial condition for / in solving the Boltzmann 
equation. As already mentioned^lhis equation only takes into account num- 
ber conserving 2 — > 2 processeal^l. 2 — > 3 processes may have a large effect 
although how significant the effect may be still remains to be determined. 
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